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We investigate the stability of a new class of BPS cosmic strings in N=l supergravity with D-terms 
recently proposed by Blanco-Pillado, Dvali and Redi. These have been conjectured to be the low 
energy manifestation of D-strings that might form from tachyon condensation after D- anti-D-brane 
annihilation in type IIB superstring theory. There are three one-parameter families of cylindrically 
symmetric one- vortex solutions to the BPS equations (tachyonic, axionic and hybrid). We find 
evidence that the zero mode in the axionic case, or s— strings, can be excited. Its evolution leads to 
the decompactification of four-dimensional spacetime at late times, with a rate that decreases with 
decreasing brane tension. 

PACS numbers: 11.27. +d, ll.25.Mj, 04.65.+e, 12.60.Jv 



Recently Blanco-Pillado, Dvali and Redi have pro- 
posed a model to describe a D-brane anti-D-brane 
unstable system after compactification to four dimen- 
sions Jl|. In the type IIB case, the tension of the 
branes appears in the four-dimensional effective theory 
as a constant Fayet Iliopoulos term which allows for 
the existence of non-singular BPS axion-dilaton strings 
generalising earlier work by in heterotic scenarios. 

It was shown that this model contains three different 
families of BPS cosmic string solutions: tachyonic or cf>- 
strings, with winding n (see also yj), axionic 1 or s-strings 
with winding m and a third type which is the special case 
where the following relation is satisfied: q\m\ = \n\. We 
are using the same notation of 1] . We call the third type 
hybrid strings. Each of the families is parametrized by a 
real number, k, defined by the relation: 

S 6 
s = 2-(\n\ -g|m|)logr-2-log|0|+K. (1) 
q q 

In the first two families the parameter k is associated to 
a zero mode that connects vortex solutions with different 
core radius and equal magnetic flux, as can be seen by 
setting: 



k = 2-{\n\ - q\m\)\ogR 

q 



(2) 



In this way R gives the scale of the core radius. In the 
third case the same parameter measures which of the 
fields, the dilaton or the tachyon, contributes more to 
compensate the Fayet-Iliopoulos term. 



1 The name axionic is somewhat misleading for the solutions dis- 
cussed here since there is no aFF coupling, in particular they do 
not share the features usually associated with axionic strings Q • 



Supergravity effects were considered in pj and, as 
expected on general grounds @, Q, the zero mode 
survives the coupling to supergravity. 

The s-strings are peculiar . The authors of [jj argued 
that they should be associated with D anti-D bound 
states that are unstable in ten dimensions, and therefore 
only exist after compactification. They also noted that 
they share some features with semilocal strings in the 
Bogomolnyi limit In the semilocal case any excita- 
tion of the zero mode leads invariably to the spread 
of the magnetic field and the eventual disappearance of 
the strings |To| . 

We have studied numerically the dynamics of this zero 
mode and we find that also in the s-string case it can 
be excited and will lead to the dissolution of the strings. 
Since in this process the field s grows without bound, 
and therefore also the compactification volume modulus, 
our result would appear to imply the decompactification 
of spacetime from four to ten dimensions at late times. 

We start by briefly reviewing the model of Q] and the 
three families of BPS strings. We then analyse the zero 
mode numerically and conclude that it can be excited. 
For comparison, we show the results of the same analysis 
on the ^-string, where a perturbation of the string leads 
to some oscillations but no runaway behaviour. 



I. THE MODEL 

The authors of Q proposed a supersymmetric abelian 
Higgs model containing a vector superfield V and N chi- 
ral superfields. Here we only need to consider one chiral 
superfield $ with charge q, that represents the tachyon. 
The model also involves an axion-dilaton superfield S 
which is coupled to the gauge multiplet in the usual way. 
Its lowest component is s + ia where the axion a is the 
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four-dimensional dual of the Ramond-Ramond two-form 
zero mode after compactification. Its scalar partner s is 
some combination of the dilaton and the volume modu- 
lus. We take the Kahler form to be K = -M 2 \og{S + S) 
and the gauge kinetic function are set to be constant, 
f(S) = I/.9 2 , where g is the gauge coupling. With this 
choice, there is no aF^F^ v term in the action. The U(l) 
symmetry is not anomalous, it is the diagonal combina- 
tion of the U(l)s from the D- anti D system. 

In component notation the bosonic sector of the la- 
grangian, after eliminating the auxiliary field from the 
vector multiplet, is: 

J2f = -\D^\ 2 - K S§ \D„S\ 2 -ig- 2 F» u F^ 

-±g 2 (t + 25K s -qW 2 ) 2 . (3) 

Ks and K s § represent the derivatives of the Kahler po- 
tential respect to the fields S and S. <f) is the lowest 
component of the chiral field field, is the £7(1) gauge 
field, and F^ — d^Ay — dyA^ the associated abelian field 
strength, 

= 0„0 - ig^V, D ll S = d ll S + i28Ar (4) 

5 is the coupling of the axion to the gauge field. In the 
case of anomalous U(l), where f(S) = S, it is called the 
Green-Schwarz parameter and is fixed by the value of 
the anomaly. Here 8 is not determined. 

It is convenient to write the bosonic part of the la- 
grangian using the rescalings: 

<t> = VUq 4> s = 5M 2 /(, s a = 25/q a 

A„ = g\/Ug K x = (gVZq)' 1 x , (5) 

With these definitions the axion a is defined modulo 
2tt, and 8 is rescaled away. After dropping the hats, the 
bosonic sector of the lagrangian reads: 

^{tg)- 2 = -\D^\ 2 -Has)' 2 (d^s) 2 - (a/ s) 2 (d^a + A^) 2 
-iF^F^-^l-s- 1 -]^ 2 ) 2 , (6) 

with D/j,(f> — dfj,cj) — iApcj), and a 2 = £/(qM 2 ). Note that, 
since a is the symmetry breaking scale in Plank units, 
ignoring supergravity and superstring corrections would 
only be a consistent approximation for a <C 1. How- 
ever, this limit is difficult to analyze numerically. We 
will present numerical results for a = 1 and argue sepa- 
rately on the effect of lowering a. 

To study straight vortices along, say, the z— direction, 
we drop the z dependence and set A z — 0. For time 
independent configurations and defining S = s + 2ia 2 a, 
and Dfj,S = d^S + 2ia 2 A ll , the energy functional can be 
written in the Bogomolnyi form: 

S{^g)- 2 = \{D X ± iD y )(j)\ 2 +±(as)~ 2 \(D x ± iD y )S\ 2 
+ h(F xv T{l-s- 1 ~W)) 2 ±F xv 
Ti[d x {<t>*D y <t>)-d y (<p*D x( />)} 
± i^a- 2 [d x {s- x D y S) - dyis-'Dj)}. (7) 



leading to a bound on the energy per unit length 

E =j d 2 xS > (£s) 2 / d 2 xF xy . (8) 

The bound is attained by the solutions of the Bogomolnyi 
equations 

(D x ± iD y )(j) = 
(D x ±W y )S =0 
F xv T (1 — s~ l — \<f>\ 2 ) -0 (9) 

We will focus on cylindrically symmetric vortices, that 
can be described by the following ansatz: 

cj) = f(r)e lnS s- 1 = h 2 

a = m0 A g = v(r)/r . (10) 

(The ansatz for the dilaton makes comparison to the 
semilocal case easier. Also we shall see that h vanishes 
in various cases, and with this choice we avoid having to 
deal with infinities). 

With this ansatz the energy density becomes: 

S = A(r) 2 + B(r) 2 /(a h) 2 +±C{rf + D(r)'/r (11) 

and the Bogomolnyi equations 

A(r) = f'-f(\n\-v)/r = 

B(r) = ti -a 2 h 3 (\m\-v)/r = (12) 

C(r) = v'/r-(l- f -h 2 ) = 

D(r) ee v- f 2 (v - \n\) - h 2 (v - |m|) = (13) 

The total energy is: 

E = 2irfdrr£'(r) > 2ir Voo (14) 

where Vqo is the asymptotic value of the gauge field profile 
function for large values of r. 

The condition A{r) = B{r) = implies the following 
relation between the tachyon and the dilaton: 

\/{ah) 2 = 2{\n\ - \m\) logr - 2 log/ + k. (15) 

Depending on the asymptotic value of the profile 
functions for large values of r, /oo, ft-oo and Woo, we can 
classify the solutions to the Bogomolnyi equations in 
three different families. Each of them is parametrized 
by the integration constant k. 

• In the first case, 

/oo = 1, = 0, Voc = n , (16) 

the tachyon acquires a non vanishing vacuum expectation 
value far from the center of the string. The magnetic flux 
of these vortices is induced by the winding of the tachyon, 
n. 

The profile function h(r) tends very slowly, (loga- 
rithmically), to zero at large r. The details about the 
asymptotics of the fields can be found in Q]. Following 



3 



to a Nielsen-Olesen string. In the case of s— strings the 
width of the string increases with decreasing a, while the 
condensate does not vary much. In fact, the main effect 
of decreasing a will be a slowing down of the dynamics. 




FIG. 1: The axion string profile functions, for m=l, n=0, 
a(r), (upper dashed line), h(r), (middle dashed line) and f(r), 
(lower solid line), with the size of the condensate, fo = 0.51. 



Blanco-Pillado et al. we call these vortices (^-strings. 

• In the second case, Fig.Q, the dilaton alone is re- 
sponsible for compensating the D-term. The function 
h(r), approaches a non vanishing constant far from the 
centre of the string, while the tachyon expectation value 
tends to zero. In this family the magnetic flux is induced 
by the winding of the axion, m. 



o h c 



i 



(17) 



These vortices have been denominated s-strings. They 
are regular thanks to the constant Fayet-Iliopoulos term 

In the previous two families the constant k has the 
same interpretation. Each value of this parameter is 
associated to a particular width of the strings. 

• For strings of the third family both the tachyon and 
the dilaton contribute to compensate the D-term. This 
happens when the axion and dilaton have the same wind- 
ing. 



f 2 + h 2 = 1 

Joe I • "oo L 



(18) 



In this case and h^ can have any value as long as 
the previous relation is satisfied (|18fl . Each particular 
/oo can be associated to a single k, which means that 
the interpretation of this parameter is different to the 
previous two cases, and cannot be related any more to 
the width of the strings, so we will not discuss it any 
further. 

As can be seen in (I13f) the derivatives of h scale as 
a 2 . In the case of the ^-strings, varying a for a fixed 
value of the integration constant k does not change much 
the width of the string, but the condensate flattens. In 
the limit when a is very small the strings are similar 



II. DISCRETIZED EQUATIONS OF MOTION 

The functions f(r), h{r) and v(r) are substituted by 
the set of quantities fk, hk and Vk which are the profile 
functions evaluated in the lattice points r k — (k + i)A, 
where A is the lattice spacing. 

We want to analyze the response of the solutions of 
the Bogomolnyi equations under perturbations. We 
must make sure that the configurations that we are going 
to perturb are stationary solutions of the equations of 
motion. This is automatic in the continuous case but 
not in an arbitrary discretization. 

Following Leese [ljj we construct a discrete version of 
the energy functional for static configurations given by: 



-Al 



±B 2 k 



a 2 (h k +i + h k ) 2 



where 

Ak = fk+i - fk - (n — v k ) 
B k 
C k 



2A 2 



(D, 



k+1 



k +J 



9) 



hk+i 

Vk+l 



hk - a (m— Vk) 

Vk 



/fc+l + fk 

2k + 1 ' 
(h k+ i + h k ) 3 



2k + 1 

J. i ; - J-2 



-A 2 (l-h 2 +1 -f 2 +1 ), 



v k - f k ( v k ~ n) - h k (v k - m). 



(20) 



The profile functions are obtained minimizing the total 
energy for static configurations: 



E = 27rJ2A 2 (k+r)g k , 



(21) 



k=0 



for which the discretized Bogomolnyi equations have to 
be satisfied: Ak = Bk = Ck = 0. 

The boundary conditions used to solve them are set 
at the points r$ and r\. At r we impose Vq = 0, and at 
n we fix the values of f\ and h\. One of these two has 
to be tuned in order to obtain the correct asymptotic 
behavior, and each value of the other one corresponds to 
a different solution within a family. 

The following discretized version of the action is nat- 
urally associated to the energy functional (|21() : 



S- 



r 2 ^ 



-2tttA- 



(fl +1 ~ fk) 



l,k=0 
I \2 



(k+i)(4-^k) 



a 2 h\ 



{2k + 1)A 4 



a ^22) 
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The superscript I labels the time slices, which are sep- 
arated by an interval r. 

Here ^ represents the density of energy associated to 
the time derivatives. 

The equations of motion can be derived from 1)22(1 by 
setting to zero the partial derivatives of the action with 
respect to the variables f l k , h l k , v l k . 



As the solutions of the Bogomolnyi equations are 
static and minimize the discretized energy functional 
they must be stationary solutions of the discretized time 
dependent equations of motion. 




The boundary conditions have been implemented us- 
ing the same method of ^3] ■ All the quantities measured 
take information from a region of radius r ca i = 4.5. The 
simulation is stopped at t max — 2(r max — r ca i). In this 
way the region from which we take data is not affected 
by the presence of the boundary. A typical value for the 
size of the lattice is r max = 21.5, but it varies depending 
on the profile. 

During the simulation we keep track of the following 
quantities: 



E(l,k ca i) 

T(l, k ca i) 
Et(1, kcal) 
W 

F(l,k ca i) 



2^A 2 ]T(fc+i)^ 

fc=0 

E(l, k ca i 



T(l,k, 



cal ) 



E(0,k cal )+T(0,k cal ) 

E(l,k ca i) 
2w(D{ cai - D l ) 



(23) 

(24) 
(25) 

(26) 
(27) 



Here fc ca/ is defined by r cal = (± + k ca i)A. 
E is the static energy contained in r < r ca i, and T 
the energy due to the time derivatives in the same 
region. Et is the total energy normalized to the initial 
value. W is a measure of the width of the string. In 
B26JI the expression in the numerator is similar to the 
energy functional, but the extra factor (fc+i) gives more 
weight to the energy far away from the core. Then, W 
increases as the energy spreads. F gives the magnetic 
flux confined in the region r < r ca ;. 

To implement the perturbation, we set a solution of 
the discrete Bogomolnyi equations in the first time slice 
I = 0, and in the second, 1 = 1, we put the same solution 
slightly deformed. We characterize the strength of the 
perturbation by the fractional change of width between 
the first, (/ = 0), and second, (/ = 1), time slices: 



FIG. 2: Response of a cj> — string (n = 1, m = 0), with 
a = 1 and condensate size hi — 0.5, to a perturbation with 
strength Wo = —0.261. The plotted lines are, from top to 
bottom, E, F, W, Et and T (except for W and Et, the rest 
of the plots have been rescaled by a factor of 1/2 to fit in the 
window). The core width, W, oscillates but is constant on 
average showing that the zero mode is not excited. 



During the simulation we have set A = 0.1 and r = 0.05, 
with r/A smaller than the Courant number l/y/2. 



A. Tachyonic Strings 

In this case to implement the perturbation we take the 
second time slice to be : 



fl = (l+p{r))f k , al = (l+p(r))a° k 
4 = (l/(^) 2 + 21og(l-Kr))r 1/2 



(29) 



Wo = (W(l = 1) - W(l = 0))/W(l = 0). 



(28) 



with the perturbation p(r) = po [1 — 3(r/r ) 2 + 2(r/r ) 3 ] 
for r < tq and zero otherwise. 

The perturbation has been chosen in order to maximize 
the fraction of energy absorbed by the zero mode. Notice 
the relation between the perturbations of the tachyon 
field and the dilaton, and the equation l|15|) that gives h 
in terms of /, in fact for small values of the parameter po 
the perturbed profile approximately satisfies the Bogo- 
molnyi equations. This perturbation gives a coordinate 
dependence to the integration constant k. 
The results shown are for a value of po such that the 
perturbation initially reduced the width of the string, 
but the same results are obtained in the opposite case. 
The parameter ro varies for different initial conditions. 
In general it has a value close to ro = 2. 

Fig.J2J shows the evolution of a 0-string, with n = 1, 
m = and a core size hi = 0.51. We show the case 
a = 1, which is also the choice made in JlJ. The 
perturbation applied has a strength Wq — —0.261, which 
corresponds to a 0.4% perturbation in the energy. We 
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B. Axionic Strings 




The relevant perturbation that excites the zero mode 
in this case is: 




ti = (l-p(r))fl 4 = (l-p(r))a° k 
h\ = (l/(h° k ) 2 + 21og(l+p(r)))- 1/2 



(30) 



FIG. 3: Response of an s— string (n = 0, m = 1), with a = 1 
and condensate size /i = 0.49, to a perturbation with strength 
Wo — —0.095. The plotted lines are, from top to bottom, E, 
F, W, Et and T (except for W and Et, the rest of the plots 
have been rescaled by a factor of 1/2 to fit in the window). 
The core width, W , after a transient, oscillates and increases 
at a constant rate, the zero mode is excited in this case. 



have plotted the observables defined in the last section 
as a function of time. Upper solid line represents E, the 
dashed line just below is the magnetic flux in r < r ca i, 
which remains almost constant. The kinetic energy, T, is 
represented by the bottom solid line. The dashed line at 
the center of the figure is the total energy. Although it 
can not be clearly seen in the plot, the data tell us that 
during the period before t = 12, a fraction of the energy 
is lost. This is the initial burst of radiation emitted 
after the perturbation. After that the system reaches a 
stationary state where all quantities oscillate except the 
magnetic flux and the total energy. 



The result of applying this perturbation with a 
strength of Wq = —0.095 to an s-string, with a = 1, 
can be seen in Fig.©. In this case the perturbation 
in energy is 2.4%. The string has windings n = and 
m = 1, and core size /i = 0.51. The functions plotted 
Timeare the same ones that appear in Fig.©. One of the 
most relevant features of these plots is that the magnetic 
flux and the total energy are decreasing with time, 
which implies that the energy is flowing out from the 
region r < r ca i, and the magnetic flux is spreading. 
At the same time the width of the string, ignoring the 
oscillatory behavior, increases at a constant rate. Before 
t = 5 oscillations are noisy. In this period the shock 
wave produced by the perturbation is still inside r < r ca i . 

Although the perturbation was chosen to reduce the 
core width, the time interval when the core is contracting 
cannot be seen clearly in the figures. The reason is 
that the contracting regime ends before the initial 
burst of radiation comes out from the observed region. 
As the system is not in a steady state yet, the data 
are difficult to interpret. We have chosen to show this 
case because the expanding regime is shown more clearly. 



The remaining solid line in the center is the string 
width. As the energy contained in the region r < r ca i 
remains constant and the width of the vortex oscillates 
only around the initial value we conclude that this 
kind of string is stable under this perturbation. The 
experiment has been repeated in a wide range of the 
parameters, p$, r$, and for different initial widths 
(parametrized by hi), and windings, but the results 
are similar to the ones presented here. Nielsen- Olesen 
vortices react in the same way to a perturbation, as was 
shown in |l(j . 

We have repeated the evolution for different values of 
a but no qualitative change has been observed. As was 
mentioned before, the smaller the value of a, the more 
similar the cf>— string is to a regular Nielsen-Olesen string, 
which is known to be stable. 




FIG. 4: Response of an n = 0, m = 2, a — 1 s-string with a 
core condensate of size /i = 0.51 to perturbations with differ- 
ent strengths. Dashed lines correspond to F, the strength is 
lowest for the top one. Solid lines represent W , the strength 
is highest for the top one. The strengths are: Wo = —0.004, 
-0.014 and -0.024. F is rescaled by a factor of 1/4. The 
figure shows how the growth rate of the core increases with 
the perturbation strength. 
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FIG. 5: Response of an n = 0, m = 1 s-string with a core 
condensate of size fi = 0.51 to a perturbation with Wo = 
0.233. The curves represent F, the magnetic flux rescaled by 
a factor of 1/4. From bottom to top the values of a are are: 
a — 0.95, 0.94 and 0.90. The figure illustrates the slowing 
down of the core expansion with decreasing a 



In Fig. (JSJ it can be seen how the rate of expansion of 
an s— string, (to — 1, n = 0), is affected by varying the 
value of a. In this case the perturbation has been chosen 
to initially increase the core size Wo = 0.233 > 0. As we 
decrease a, keeping the perturbation strength fixed, the 
rate of expansion of the string decreases. This can be 
understood from equation (|22l) . The energy associated 
to the field h scales as the inverse of a 2 . Deviations from 
the solution to the Bogomolnyi equations cost more 
energy for smaller values of a, thus for a fixed perturba- 
tion strength the evolution rates should decrease with 
a. The values of alpha are: a = 0.95, 0.94, 0.90. 

The precision of the technique used here does not 
allow to obtain reliable data for values of a lower than 
0.7, where already the evolution is so slow that it can 
hardly be appreciated during the time of the simulation. 
However, note that in these simulations time is measured 
in units of the inverse of the Higgs mass, thus even for 
lower values of a decompactification is still possible on 
cosmological time scales. 



Fig. |0| shows the effect of applying perturbations of 
different strengths toan = 0,m = 2 vortex. In this case 
the vortex has also a condensate size of /i = 0.51. The 
strengths applied are: Wo = -0.004, Wo = -0.014 and 
W = -0.024. 

Notice that the bigger the strength of the perturba- 
tion, the larger the fraction of magnetic flux lost through 
the boundary r = r ca i . The rate of growth of the radius 
also increases with the strength of the perturbation. 
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